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We consider a SBS Josephson junction the superconducting electrodes S of which are in contact 
with normal metal reservoirs (B means a barrier). For temperatures near T c we calculate an effective 
critical current /* and the resistance of the system at the currents I < I* and I » I*. It is found 
that the charge imbalance, which arises due to injection of quasiparticles from the N reservoirs into 
the S wire, affects essentially the characteristics of the structure. The effective critical current /* 
is always larger than the critical current I c in the absence of the normal reservoirs and increases 
with decreasing the ratio of the length of the S wire 2L to the charge imbalance relaxation length 
Iq. It is shown that a series of peaks arises on the I — V characteristics due to excitation of the 
Carlson-Goldman collective modes. We find the position of Shapiro steps which deviates from that 
given by the Josephson relation. 

PACS numbers: 74.20.De, 74.25. Ha, 74.25.-q 



I. INTRODUCTION 

In the theory of the Josephson effect in weak links [l( it is assumed that the superconducting electrodes are in 
equilibrium. This means, in particular, that a gauge-invariant potential /i 

li=\d X /dt + eV, (1) 
is equal to zero [2] (we set h = 1). The Josephson relation follows immediately from this condition 



2e(V 2 -V 1 )=d<p/dt, (2) 

where (V2 — Vi) is the voltage drop across the Josephson junction and <p = xi ~ Xi is the phase difference. Most 
works on studies of the Josephson effects were carried out under equilibrium conditions so that the potential \i is zero 
in the superconducting electrodes S\fi, and Eq. ^ is satisfied The potential /i is related to a so called charge 

imbalance arising due to a different population of the electronlike and holelike branches of the quasiparticle spectrum 
of a superconductor I5|. The conditions under which this charge imbalance may arise were studied experimentally 
@, 0, i, 1 E3, El El, HI and theoretically 0, El El E3, El, El ■ 

In the last decade a great progress in preparation of superconductor/normal metal (S/N) nanostructures has been 
achieved. Varies properties of these structures were studied experimentally. One can mention the study of transport 
[l2l , I20L l2ll [221 . [23l . Izij and thermoelectric [H [27j properties, the measurements of the density-of-states [l^ etc. 
Recently a transition of a thin and short superconducting wire into a resistive state caused by a current I was 
investigated [28|. In particular, an increase of the critical current I c has been observed when an external magnetic 
field H is applied [28[ . A possible reason for the observed increase of I c may be a pol arization of magnetic impurities 
by the magnetic field H [28|, [2^] . Another mechanism has been suggested in Ref . [30|| . It has been proposed that an 
increase of the critical currents characterizing phase-slip centers in a thin S wire may be due to a shortening of the 
charge imbalance relaxation length Iq under the influence of H. However, it remains unclear how this mechanism can 
work in a pure superconducting state when there is no charge imbalance. As is known (for review see (l9l l32l |33| ) . 
the charge imbalance arises only under nonequilibrium conditions. For example, this mechanism may be realized in 
a superconductor-normal metal {SN) nanostructure with a bias current, where the c harg e imbalance is built up due 
to injection of quasiparticles into the S wire from the N reservoirs 0, EH, El, E3, EM Ell ■ 
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FIG. 1: Schematic view of the system under consideration. 



In the present, paper we investigate the Josephson effect in a NSBSN structure under nonequilibrium conditions, 
in the presence of a bias current / flowing through SN interfaces. In this case the potential fi arises due to injection 
of quasiparticles into the superconductors S from the normal reservoirs N {B is a barrier of an arbitrary type). We 
show that if the length of the S wire 2L is less than or comparable with the charge imbalance relaxation length Iq , the 
critical current increases and may significantly exceed the critical current I c in the absence of the normal reservoirs. 
The resistance of the system also will be calculated for the currents I < I* and I I*, where I* is an effective critical 
current. We find the position of Shapiro steps and show that peaks associated with the excitation of collective modes 
of the Carlson-Goldman type arise on the I — V characteristics. 

Note that in recent publications [H, [35| a NSN system without the Josephson junction was considered and the 
I — V curve of the system was calculated numerically. 



II. MODEL AND BASIC EQUATIONS 



We consider the structure shown in Fig.l. We assum that the barrier is located in the middle of the S wire, but 
the results remains qualitatively unchanged for a system with the barrier located at the S/N interface. The system 
consists of a superconducting wire (S wire) connecting two normal reservoirs N. In the middle of the S wire there 
is a barrier which provides the Josephson coupling. The transparency of the SN interface is assumed to be high. 
The analysis can be easily generalized to the case of arbitrary NS interface transparencies. We also assume that the 
temperature T is close to the critical one T c and therefore the inequality 



A«T (3) 

is satisfied (A is the order parameter in the S wire). In this case the effects of the branch imbalance are most 
significant. The lengths L and Iq are assumed large in comparison with the Ginsburg-Landau correlation length £gl, 
i.e., 



£,gl « l.2^D/T c (T c /A) « {L, l Q }, (4) 

where D is the diffusion coefficient in the S wire. The charge imbalance relaxation length Iq is determined by inelastic 
scattering processes and may be rather long ( fl5l. [lil [l7l. pH Il9j ) . This assumption allows us to consider the order 
parameter constant in the major part of the S wire. The same limit was considered in Refs. (l7t HH where the 
resistance of the superconductor in a SN system was calculated. The relation between L and Iq may be arbitrary. 

Our aim is to obtain a relationship between the current I and the voltage Vl = V(x) at x = L (the voltage difference 
between the N reservoirs is Vnn — — 2Vt). We restrict ourselves with voltages Vnn small compared to the energy 
gap A : eV/vw ^ A. In this limit the distribution function (longitudinal in terms of Ref. [Hj]), which determines 
the order parameter A, is close to the equilibrium one. Another distribution function denoted by fi in Ref. (3l| 
(transverse in terms of Ref. [Ill) was found in Ref. [H, To be more exact, the so called anomalous Green's 
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function, = g R fsfi — fifsg A , was found in Refs. [13, [H, Oil , where g R ^ are the retarded (advanced) Green's 
functions. Using the functions g^ and g R ^ A \ one can obtain the expression for the current / in the 5* wire, which in 
the main approximation in the parameters A/T and V/T is equal to 

1 = Sa{E + W~e Ql (5) 

where S is the cross section area of the S wire, a is the conductivity of the S wire in the normal state. We assume that 
the cross section area is small compared to {£,ql, A^} so that all the vectors {I, E, Q} have only the x— component 
{/, E,Q} and depend on x. The electric field E = —dV(x)/dx (we drop the vector potential A because we consider 
only the longitudinal electric field choosing a corresponding gauge). The condensate momentum Q is defined as 



Q = (l/2)dx/dx~2nA/^ (6) 

where $o = hc/2e is the magnetic flux quantum; the vector potential A can be dropped because we do not consider 
the action of magnetic field. The momentum Q obeys the equation 



dQ/dt=eE + d/J l /dx, (7) 
The spatial and temporal variation of the gauge-invariant potential fJ.(x,t) is described by the equation (l9l. 1371] 



(d/dt + 7 )/i = v 2 CG dQ/dx, (8) 

where 7 is a quantity which determines the charge imbalance relaxation rate. It is related to inelastic scattering 
processes (7 = 1/t c , where r e is the inelastic relaxation time), the condensate momentum (7 ~ DQ 2 A/T) in the 

. The velocity vcg — v2DA is the velocity of 
. In the stationary case Eq. §8§ can be written 
in the form 



presence of condensate flow or the gap anisotropy [1 
propagation of the Carlson-Goldman collective mode 11 
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Z Q V = (e/a)d 3s /dx = d 2 ^/dx 2 , (9) 



where js = a(nA 2 /2Te)Q is the density of the supercurrent; Iq = y/ ATD/ (njA) is the penetration depth of the 
electric field (or the charge relaxation length). Eq.® describes the conversion of the quasiparticle and superconducting 
currents. One can see that the potential /i arises if the divergence of the supercurrent (or quasiparticle current) differs 
from zero. 

The current in the S wire can also be written as the current through the Josephson junction. In the main approxi- 
mation it is equal to (see Appendix) 



I = — -^+/ c sin^, (10) 
Kb 

where 2Vq = 2V(0+) = — 2V(0— ) is the voltage drop across the Josephson junction (Vq is negative), Rb is the barrier 
resistance, I c is the Josephson critical current, and ip = x(0+) — x(O-) is the phase difference. For simplicity we 
drop here the displacement current CdVo/dt assuming that the parameter (Ci?B)e(/ c -Rs) is small (the main results 
concerning the critical current and the resistances of the system do not depend on the presence of the displacement 
current). At V(L) <§; A the critical current equals I c = irA 2 /ATeRs- The first term in Eq. tflT))) is the quasiparticle 
current and the second term is the supercurrent. 

Eqs.([5] fTU|) describe the system. These equations should be complemented by boundary conditions. Since we consider 
the voltages smaller than A, the distribution function which determines the supercurrent, is close to the equilibrium 
one: /; = tanh(e/2T). This implies the conservation of the quasiparticle (correspondingly superconducting) currents 
at the Josephson junction, i.e., 



ScrE^O-i-) = — — — = I — I c s'mip. 
Rb 



(11) 
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Another boundary condition relates the electric field at the edge of the S wire E(L) to the electric field in the N 
region Epf. As is shown in Refs. (TtI Il8j ■ the electric field E(x) experiences a jump at the SN interface 

[E] SN = E N - E(L) = rE(L) (12) 

where r w 0.7((A/T)IJCgl) 1/4 « 1.17[((T , c -T)/T c ) 2 t £ /t] 1 / 4 ; t £ ,t are the inelastic and elastic scattering times. This 
jump is small if the temperature T is close to T c and therefore the electric field (accordingly the quasiparticle current) 
is continuous at the SN interface. For a finite value of r and equal conductivities in the S and N regions we get from 
Eq.® 

E{L) = E N {1 - (13) 
1 + r 

Here the second term on the right is small near T c . Solving Eqs. (|5l7l8ll0p with the boundary conditions (|1HI13[) , we 
can calculate, in particular, the effective critical current I* and the resistance of the system R. 

III. RESISTANCE AND JOSEPHSON CRITICAL CURRENT. 



Eqs. (|5l7l8ll0p are nonlinear equations in partial differentials in which all functions depend on time t and coordinate 
x. Therefore the solution to these equations can not be obtained in a general case. We consider limiting cases of small 
and large currents I. 

a) Stationary case; I < I*. 

If the current does not exceed a critical value 7* (the effective critical current I* will be found below), the phase 
difference ip and other quantities do not depend on time. In particular, dQ/dt = and the electric field is 

eE(x) = -dfi/dx (14) 
The momentum of the condensate is found from Eq. §5§ 

vrA 2 

— Q = eIR n + d[i/dx (15) 

where R n = (So~)~ l is the resistance of the S wire in the normal state per unit length. The potential /i is described 
by Eq.((5]). At small currents when the condensate flow does not contribute to the relaxation of the charge imbalance, 
the length l Q equals: l Q = ^4(DT e )T/nA [II []]], [HI . The solution of Eq.® is 

fi(x) = Acosh(x/l Q ) + Bsmh(x/l Q ) (16) 

The electric field and the potential V(x) equal 

eE(x) = -[Asmh(x/l Q ) + B cosh(x/l Q )]/l Q , eV(x) = Acosh(x/l Q ) + Bsmh(x/l Q ) + K (17) 

where V* is an integration constant. 

First we neglect the second term on the right in Eq. (TT3"]) . From Eqs. (|llll3[) we find the coefficient B and A 

B = —cRq (I — I c sin ip) , A= -(eR Q I + BC)/S (18) 

where Rq = Iq/(S(t), C = cosh^ and S = sinh^ with 9 — L/Iq. 
From Eq. (fT7|) we obtain 

eV L = eV Q + A(C -l)+BS (19) 
Writing Vq in terms of B (see Eqs. (|llll8p ), we can represent Vl in the form 
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eV L = S(~ + ^i) - IRq——^— (20) 
In the stationary case one has eVb = /i(0) = A. Combining Eqs. (|llll8[) . we obtain from this formula 

^+^)=^4 + § )sin ^ (2i) 

where b = Rb/Rq. Eq. ([2T|) determines the effective critical current 7* 

r / 6S + 2C _ 6tanh(0/2) + l + tanh 2 (g/2) 

c/ c bS + 2(C-l) tanh(6»/2)(6 + 2tanh(6»/2)) 1 j 

It is seen that the effective critical current I* is always larger than the critical current I c in the absence of the charge 
imbalance. In the limit of large b = Rb/Rq or 9 — L/Iq the effective critical current coincides with I c . Therefore in 
a Josephson junction with a weak coupling between superconductors (b 3> 1), the effects of charge imbalance are not 
important. For a short S wire, we obtain 

r/j- b9 + 2 (23) 

1Jlc ~ 9{b + 29) {26) 

One can see that I* diverges at 6 — > 0. This fact is in agreement with the results of Ref.[4(|, where a very short 
NSS' system was considered, and it was concluded that the stationary state exists at any currents I. However our 
consideration is valid only for not too small 8 (9 = L/Iq ^ £,gl/Iq)- The dependence of I*/I c on 9 is shown in Fig. 2 
for different b. 

What is the reason for the increase of the effective critical current from the physical point of view? The critical 
current is defined as a maximum current I ma x at which the stationary state is possible: dtp/dt = 0. The current I max 
in an ordinary equilibrium Josephson junction coincides with the critical current I c because the first term in Eq. (|10j) 
in the stationary state is zero: -2Vq = (h/2e)dip/dt = 0. As we noted in the Introduction, this equation follows from 
the fact that the potential fi(x) is zero. The same situation takes place near the S/B/S interfaces in a long (L >> Iq) 
S wire. The charge imbalance, which determines fJ.(x) arises only in the vicinity of S/N interface and decays on the 
scale of order Iq. Therefore the effective critical current is the maximum current at which the stationary state with 
nonzero potentials V and /i is possible. Here we find this current for the case of a short system: L << Iq. Then, the 
electric field E(x) and potential V{x) can be written as 

eE(x) S -[A{x/l Q ) + B]/l Q ,eV{x) = A + B(x/l Q ) (24) 



The first term in the square brackets of Eq. (l24[) describes the conversion of the quasiparticle current into the 
condensate current. In the normal state it is zero because Iq — > oo. At x = L the charge is transferred by 
quasiparticles (we ignore the jump in the electric field setting r — 0): SaE(L) = I. At x = the quasiparticle current 
is SaE(0 + ) = I — I c smif. Therefore the coefficient A related to the conversion of the quasiparticle current into the 
condensate one is equal to 

- A(L/l Q ) = eI c R Q sin <p. (25) 
On the other hand, at currents I large in comparison with I c one has 

I ^ -2Vo/R B = -2A/eR B - (26) 
Thus, from Eqs. (|25H26p . we find the maximum current for the stationary state 

r c ^I c (2l Q /bL) (27) 
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FIG. 2: The normalized effective critical current as a function of 9 = L/Iq for b = Rb/Rq = 0.2 (curve 1), b = 2 (curve 2), 
and b = 5 (curve 3) 

This formula is valid if the ratio (Iq/bL) is large. If the current exceeds this value, the stationary state is not 
possible. One can say that the length of the S wire is too short to provide the conversion of the quasiparticle current 
into the condensate current I c simp. 

If we take into account the second term in the boundary condition (|13p , we obtain for I* 



2 + 6 tanh f 



tanh 9[b + 2 tanh(0/2) +2r/(S(l + r))] 
The resistance of the system R = Vnn/I = |2Vl|/7 at I < I* can be easily found with the aid of Eqs. (|20|21|) 



(28) 



Rn — 2R, 



2 tanh 6 + b 
'2 + b tanh 6 



(29) 



In limiting cases we find 



. (28 + b)/(2 + b8), 0«1 
Ro - 2R Q ^ g >> 1 



(30) 



The resistance Rq is equal to Rb at 9 —> and to 2Rq at 9 — > oo. Thus, the resistance of a short system is a 
combination of the barrier resistance Rb and Rq (the resistance of the S wire near the SW interface). The resistance 
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FIG. 3: The resistance of the system at currents less than the effective critical current. as a function of 9 for b = 0.2 (curve 1), 
b — 2 (curve 2), and b = 5 (curve 3) 

of a long system equals the resistance of the S wire near the SN interface on the scale Iq. The dependence of the 
resistance Rq on 8 for different b is shown in Fig. 3. 

With account for a hnite jump of the electric field at the SN interface (the second term on the right in Eq. (fTS"!) ). 
we obtain for Rq 



Rq 



1 



2 2tanh6> + 6 
~+r Q 2 + 6tanh0 



(31) 



Consider now the case of large currents I when the phase difference cp(t) is increasing in time (<p(t) ~ t) and its 
oscillating part is small. 

a) Quasistationary case; I I*. 

In this case the condensate momentum is almost time- independent so that dQ/dt = 0. The potentials //, V and the 
electric field E are described by Eqs. (I16I17[) . but the formula for the coefficient B is changed: B = —bRqI. From 
Eq. ([2U)) we find the voltage difference between the N reservoirs 



V NN = -2V L = I[R B + AR Q < ^-] (32) 



and the resistance Roo at large currents 



R x = [R B + 4R Q tanh(6>/2)] 



(33) 
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In the case of a long S wire (9 ^> 1) we obtain: i?^ = Rb + 4-Rq; that is, the resistance of system is the sum of 
the barrier resistance and the resistance of the regions of the S wire where the electric field penetrates (close to the 
SN interface and to the barrier). In the case of a short S wire (6 <C 1) the resistance is: i?oo = Rb + 28Rq, that is, 
the contribution of the superconducting regions to the resistance decreases. 

With account for the second term in Ea. (fl3|) the resistance i?oo acquires the form 

Roc = [Rb +*Rq(1- r + ) tanh(0/2)] (34) 
That is, the contribution of the S region near the SN interface to the resistance decreases. 



IV. THE I-V CHARACTERISTICS AND SHAPIRO STEPS 

In this Section we analyze the form of the current-voltage (I-V) characteristics and calculate the voltage Vsh which 
determines the position of the first Shapiro step. At a finite value of the ratio Rq/Rb the voltage Vsh differs from 
that given by Eq.([2]). We restrict ourselves with the limit of large currents / employing an expansion in the parameter 
I c /I. In the considered non-stationary case all the quantities depend on time, and in the lowest approximation in 
the parameter I c /I this dependence is determined by terms of the type exp(iQjt), where f2j is the frequency of 
the Josephson oscillations. Eqs. ()8l9|16|17|) can be easily generalized for the non-stationary case. The potential /i is 
described again by the equation 

d 2 fi/dx 2 = k n fi (35) 
where k n = (ifl + -7) (ill + Ha)/vq G , SIa = 7rA 2 /2T. The solution for this equation is 

^i(x) = An cosh(fcfix) + Bq sinh(fcnx) (36) 
The electric field and potential have the form 

e£(s) = - " A dn/dx, eV(x)= " A M + eK, (37) 

til + Ha lH + "A 

In deriving Eq. (|37|) we assume that the current / does not depend on time (no external ac signal). From the 
boundary conditions (|11I13|) wc find 

Bn = -eRu — ± " A (I - I c sin tp) , An = -(eR Q I(l + r)- 1 + B n C n )/Sn (38) 

i^A 

where Rq = (fc^cr) -1 , Cq = cosh(fcoL), Sn — sinh(fcnL). The electric potential at the N reservoir is 

til + Ha 2Rq 

Here we took into account the finite value of r. The equation for the phase difference tp is obtained from the 
definition of (Xq — (l/4)dtp/dt + eVo. It has the form 

L ( % )a + [Rb + 2Rn ^^ (sta , )n . (Bs + 2fl0( £|_i + _L_]7 (40) 



We rearrange the terms in Eq. (|40p and rewrite it in a dimcnsionless form 

/ dtp \ r . r- 2 . Co — 1 r , . 2 Ccj iu + loa , , . » , , x 
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where j = I/I c , = ilj / (2eI c R B ), r = t(2eI c R B ) are the dimensionless current, Josephson frequency and time, 
respectively; = b^, b L = W, 6^ = 9 ^(iu + cj A )(icj + w 7 ), # = (L/v C g)(?I c Rb) = l.ll(A/T) 3 / 2 y/T/E^, 
£ T/l = D/L 2 , uj a = n A /(2eI c R B ) = 1, w 7 = 1 /{2eI c R B ). 

The solution is sought in the form <^(r) = <ySo( T ) + ^i(t) + •••) where f n (r) ~ j _n , for n > 1. In zero order 
approximation we obtain from Eq.(|4ip 

Vq =uj Q t, u JQ =j[l + h-9-ll + r )] (42) 
o Sq (l + r)5Q 

The correction ifi(r) is equal to 

<Pi(t) = o;y 1 [i?e£(ijj) cos(wt) — ImC(u) sin(wT)] (43) 

where = [1 + (2/b LU )(C LJ /S^)(iuj + u a )/wa]. 

Knowing </?i(t), one can find a correction to the I-V curve due to the Josephson oscillations. From Eq. (|39p we find 
for the dimensionless voltage drop < um >= — < 2Vl > /(2I C R B ) 

< v NN >= j[l + \{%^ + 2(1 + r )S Q )] + l tanh ^ 2 )] < sin ^ > ( 44 ) 

Here < sin ip >=< ipi (t) cos(o;r) >= (2loj^ 1 Re C(lo). The first term is the Ohm's law at large currents. The 
second term is a contribution from the Josephson oscillating current. Therefore the correction to the Ohm's law is 
equal to 

<5v NN >=-—!— [l + ~tanh(0/2)]i2e£(u;) (45) 
2wjo b 

In Fig. 4 we plot the dependence of this correction < Svnn > on the normalized frequency ujjo, which is proportional 
to the normalized current j (see Eq. (|42|) ) for two values of damping 7. In plotting Fig. 4 the following values of param- 
eters are taken: D = 50cm 2 /s, T c — 1.3K, A/T = 0.3, b — 10, L — 0.5mkm, r e = 2 • 10~ 10 s. For these values we have: 
^ = (2eI c R B r e )- 1 = (2/TT)(T e T)- 1 (T/A) 2 Si 0.27; 8 = (it /2V2)(A 2 /T)L/VDA Si 1.11(A/Tf/ 2 (T /E Th ) l l 2 Si 0.46, 
l Q Si 2mkm, 9 = L/l Q ^ 0.25, where E Th = D/L 2 . 

It is seen that this correction has a series of peaks. These peaks are related to excitation of a collective mode of 
the Carlson- Goldman type in the system [3(| H3, [H, [39[ . The excitation occurs if the frequency of the Josephson 
oscillations f2j = ujjo(2eL c R B ) coincides with the frequency of the Carlson-Goldman modes £Icg = vcg/L. In this 
case the quantity ReC(uj) ~ tanh _1 (fcnL) has a peak. The possibility to observe such modes in another Josephson 
system was studied in Ref.[42j]. 

Shapiro steps on the I-V characteristics arise if in addition to dc current an ac current flows in the system. In the 
presence of the ac current, I~(t) = 1~ sui(Q ex t), a new term appears on the right in Eq. (|40p , which is proportional to 
1^ (t) . The position of the first Shapiro step is determined by the equation [1, 13| 

fij = (2eI c R B )cuj = Q ex (46) 

where f2j is the frequency of the Josephson oscillations. We obtain the frequency of the Josephson oscillations, uj, 
from Eq. (141|) . At large currents j in the main approximation, loj is given by Eq. (|42p . On the other hand, the 
normalized voltage 2 < vnn > corresponding to the same current j can be easily obtained from Eq. 



2 < vnn >= j[l + £(tanh(0/2) + 1 )] (47) 

Therefore the deviation of 2 < vnn > from the value given by the Josephson relation is given by the formula 

= 2 < v NN > -uj ex = 2 tanh(6>/2) - r/S 

Sh ~~ u ex (l + r)6 + 2[tanh(0/2)+r/((l + r)S)] [ ' 
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FIG. 4: Correction to the I-V characteristics for o; 7 = 0.1 (curve 1) and u; 7 = 0.5 (curve 2). 

with uj ex = lojq. One can see that this deviation can be both negative and positive depending on the parameters 9 
and r. In the limit of large b the deviation Svsh = 0, i.e. the Josephson relation is fulfilled. At an arbitrary b, the 
parameter 5vsh depends on the quantities 9 and r. Thus, by measuring the deviation from the Josephson relation 
and the resistances Rq and Roo, one can determine the ratio b = Rb/Rq, the parameter r characterizing the jump of 
the electric field at the SN interface, and the charge relaxation length Iq — L/9. 

V. CONCLUSIONS 

On the basis of a simple model we have analyzed the Josephson effects in a NSBSN nanostructure at temperatures 
close to T c . It turns out that the charge imbalance taking place in this system essentially changes the characteristics 
of the Josephson junction. If the barrier resistance Rb is not too large in comparison with the resistance of the S wire 
Rq on the scale of the charge imbalance relaxation length Iq, the Josephson critical current increases (see Eq. (|23]) ) 
and the positions of Shapiro steps deviate from their position in equilibrium Josephson junctions (see Eq. (l48p ). We 
have also calculated the resistance of the system for small (Rq for I < I*) and large (i?oo for / >> /*) currents. The 
ratio of these resistances equals 

Roo _ (b + 4tanh(6>/2))(2 + frtanhfl) 

~R^~ ~ 2(6 + 2tanh6») ^ ' 

Thus, by measuring the ratio Rqo/Rq and the deviation Svsh, one can determine the coefficient b and the the charge 
imbalance relaxation length Iq = L/9. 

We have also shown that there is a series of peaks on the I — V characteristics of the system related to excitation 
of the Carlson-Goldman collective mode. Note that even if the barrier resistance Rb is ten times larger than the 
resistance Rq, the peaks in Fig. 4 caused by collective mode excitations are clearly visible. 

Note that we have adopted several assumptions. One of them is that the voltage Vl should be smaller than the 
energy gap A. If the voltage is higher than A, new effects may arise in the system. In particular, the critical current 
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I c may change sign [43j, |4J, |45j, l46| . It would be of interest to measure the I — V characteristics of the considered 
system in a wide range of the applied voltages. 
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VII. APPENDIX 

Here we derive the boundary condition (fTTj) and the expression for the current (fit?)) . The latter formula in our case 
is not obvious because the potential /i differs from zero and the Josephson relation J2]) is not valid. Consider the 
boundary condition at the barrier B for 4*4 matrix quasiclassical Green's function g [40l l4lj | 

a(gdg/dx)^(2SR B )- 1 [g 0+ ^} (50) 

where <?o± = #(±0) are the Green's functions on the right and left from the barrier B. The (11) and (22) elements 
of the g matrix are the retarded and advanced Green's functions g R ( A > and the (12) element is the Keldysh function 
g K . If the voltage in the system V is small in comparison with A, the distribution function /; which determines the 
energy gap A is close to the equilibrium one (/; ~ tanh(e/2T)). 
The current in the system is given by 

I=\SCJ j deTr{T 3 [g R ie,t-,r)dg K (e,t-,r)/dx + g K (e,t ]r )dg A (e,t-,r)/dx}} (51) 

If we calculate the current with the help of Eqs. (|5QI51[) . we obtain Eq.©. Let us calculate the current through 
the Josephson junction using the right hand side of Eq. (|50[) . This current consists of three terms: the quasiparticle 
current I qp , the Josephson current Ij and the interference current Ii n t [H, HI- The Josephson current is 

7 -> = |£ E Tr iMUf- /-/+]} (52) 



where f± = /(±0) = (t 2 cos(y>/2) ± fi sin(<^/2))A/ \J{uj\ + A 2 ) are the Gor'kov's quasiclassical Green's functions on 
the right (left) from the barrier. Calculating the sum, we obtain 

I j = I c sin ip (53) 
with I c — ir A 2 /ATeRs- The interference current is given by 



hnt ~ 16R B 



J de(f3V )(f? + f?))(f R + f A )co S h-\eP) 



(54) 



In the main approximation Ii nt ~ Vo A/ (AT Rb)S cos iphi(A/T), where T is a damping in the spectrum of the 
superconductor. We see that the interference current is small in comparison with the Josephson current if the voltage 
Vq is smaller than the value A/ln(A/r). The quasiparticle current I qp consists of two parts: I qp — Iqp + 5I qp , where 
I qp ^ is determined by the formula 

C = / de(f3d X+ /dt)(g* - gi))(g R - g A ) cosh- 2 ( £ /3) = -R B \d X+ /dt) (55) 
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The term f3dx+/dtcosh~ 2 (ef3) is obtained from the equilibrium distribution function after the transformation of 
the Green's functions [l||: g — ► UgU + , where U = exp(if3x/2). The correction 5I qp is determined by the response of 
the system to a perturbation of the potential fj, [l^] . It is equal to 

1 f . tanh((e - Q/2)/3) - tanh((e + n/2)/3) 2^ 

= "2^ 7 de n /i0 = 7r b (56) 

where /Uo = (l/2)dxo/dt + cVq. Thus, for the quasiparticle current we obtain the first term in Eq. lfTU]) . 

The boundary condition ([50]) may be written for the retarded (advanced) Green's functions g RtyA \ Since the dis- 
tribution function /; approximately coincide with the equilibrium function, one can easily obtain from this boundary 
condition the equation of continuity of the condensate current at the Josephson junction. Therefore, the quasiparticle 
current also is continuous at the SBS junction. 
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